1
理论基础与适定性
MATH007Lesson 5
00:00
在使用四阶龙格-库塔法或亚当斯-穆尔顿公式等数值求解器之前,我们必须提出一个根本性问题: 解真的存在吗?是否稳定? 初值问题(IVP)的理论基础提供了数学上的“绿灯”,确保我们的离散化方法收敛到有意义的物理现实,而非数值噪声。

基石:利普希茨连续性

为了控制误差的传播,我们需要一个函数 $f(t, y)$,它不会出现过大的跳跃。这一性质由 利普希茨条件

定义 5.1:利普希茨条件

若在集合 $D \subset \mathbb{R}^2$ 上存在常数 $L > 0$,使得函数 $f(t, y)$ 关于变量 $y$ 满足利普希茨条件,则有:

$$|f(t, y_1) - f(t, y_2)| \le L|y_1 - y_2|$$

对所有 $(t, y_1), (t, y_2) \in D$ 成立。该常数 $L$ 是函数垂直变化的“速度限制”。

示例 1:分析利普希茨常数

考虑函数 $f(t, y) = t|y|$ 在区域 $D = \{(t, y) \mid 1 \le t \le 2, -3 \le y \le 4\}$ 上的情况。根据中值定理(或绝对值的性质):

$|f(t, y_1) - f(t, y_2)| = |t|y_1| - t|y_2|| = |t| \cdot ||y_1| - |y_2|| \le |t| \cdot |y_1 - y_2|$。

由于在该区域内 $t$ 的最大值为 2,因此利普希茨常数为 $L=2$。

区域几何完整性

我们无法在充满“洞”的区域中求解初值问题。我们需要 凸性

定义 5.2:凸集

若对于任意两点 $(t_1, y_1)$ 和 $(t_2, y_2)$,其连线段满足:

$$((1 - \lambda)t_1 + \lambda t_2, (1 - \lambda)y_1 + \lambda y_2)$$

其中 $\lambda \in [0, 1]$,且该线段也包含在 $D$ 内。这确保了求解路径不会“离开”有效计算区域。

存在性与唯一性定理

当这些条件同时满足时,我们引用 定理 5.4:若 $f$ 在凸集 $D$ 上连续并满足利普希茨条件,则初值问题 $y' = f(t, y), y(a) = \alpha$ 具有 唯一 解 $y(t)$。这解释了为何像欧拉法($w_{i+1} = w_i + h f(t_i, w_i)$)这样简单的算法,以及预测-校正逻辑这样复杂的算法都是合理的:

$WC = w_{i-1} + \frac{h}{24}[9f(t_i, WP) + 19f(t_{i-1}, w_{i-1}) - 5f(t_{i-2}, w_{i-2}) + f(t_{i-3}, w_{i-3})]$。

🎯 核心原则:适定性
一个问题如果 适定 则称为适定的,前提是存在唯一解且解连续依赖于初始数据。若利普希茨常数 $L$ 极大,问题将变为“刚性”。在刚性方程中,瞬态部分迅速衰减,但其导数(量级为 $c^n e^{-ct}$)却不会快速衰减,因此需要 算法 5.8:梯形法结合牛顿迭代 以维持稳定性。